THE PROUHET-TARRY-ESCOTT PROBLEM FOR GAUSSIAN 

INTEGERS 



TIMOTHY CALEY 

Abstract. Given natural numbers n and k, with n > k, the Prouhet-Tarry- 
Escott (pte) problem asks for distinct subsets of Z, say X = {xi, . . . ,x n } and 
Y = {yi, . . . , y n }, such that 

4 + •■■ + <= y\ + ... + y l n 

for i = 1, . . . , k. Many partial solutions to this problem were found in the late 
19th century and early 20th century. 

When n = k — 1, we call a solution X = n -i Y ideal. This is considered to 
be the most interesting case. Ideal solutions have been found using elementary 
methods, elliptic curves, and computational techniques. In 2007, Alpers and 
Tijdeman gave examples of solutions to the PTE problem over the Gaussian 
integers. This paper extends the framework of the problem to this setting. We 
prove generalizations of results from the literature, and use this information 
along with computational techniques to find ideal solutions to the PTE problem 
in the Gaussian integers. 



1. Introduction 

The Prouhet-Tarry-Escott problem, or PTE problem for short, is a classical num- 
ber theoretic problem: given natural numbers n and k 7 with k < n, find two distinct 
subsets of Z, say X — {x\, . . . , x n } and Y = {yi, . . . , such that 

n n 

(1.1) =X>" for j = 1,2,..., fc. 

i=l i=l 

A solution is written X =k Y, and n is its size and k is its degree. The maximal 
nontrivial case of the PTE problem occurs when k — n — 1. A solution in this case, 
say X = n -i Y , is called ideal. 

For example, {0, 3, 5, 11, 13, 16} =5 {1, 1, 8, 8, 15, 15} is an ideal PTE solution of 
size 6 and degree 5 since 

+ 3 + 5 + 11 + 13 + 16 = 48 =1 + 1 + 8 + 8 + 15 +15 

2 + 3 2 + 5 2 + ll 2 + 13 2 + 16 2 = 580 = l 2 + l 2 + 8 2 + 8 2 + 15 2 + 15 2 

3 + 3 3 + 5 3 + ll 3 + 13 3 + 16 3 = 7776 = l 3 + l 3 + 8 3 + 8 3 + 15 3 + 15 3 

4 + 3 4 + 5 4 + ll 4 + 13 4 + 16 4 = 109444 = l 4 + l 4 + 8 4 + 8 4 + 15 4 + 15 4 

5 + 3 5 + 5 5 + ll 5 + 13 5 + 16 5 = 1584288 = l 5 + l 5 + 8 5 + 8 5 + 15 5 + 15 5 . 
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Similarly, for a, b, c, d £ Z, 

{a + b + d, a + c + d, b + c + d, d} ~2 {a + d,b + d,c + d,a + b + c + d}, 

is a family of PTE solutions of size 4 and degree 2 due to Goldbach. In fact, this 
example was also found by Euler for the case when d = 0. Many other elementary 
solutions can be found in |13j . 

The PTE problem is interesting because it is an old problem with both algebraic 
and analytic aspects, and also has connections to other problems. Ideal solutions 
are especially interesting because of their connection to problems in theoretical 
computer science [2] and combinatorics |16j . a conjecture of Erdos and Szekeres 
fl9 l 112 ) . [3] Chapter 13], as well as the "Easier" Waring problem, which we discuss 
below. 

Given an integer fc, the "Easier" Waring problem asks for the smallest n, denoted 
v(k), such that for all integers to, there exist integers x%, . . . , x n such that 

±x\ ± . . . ± X k n = TO, 

for any choices of signs. This problem was posed by E. M. Wright as a weakening 
of the usual Waring's problem, which allows only addition. Note that v(k) is 
conjectured to be O(k). For arbitrary fc, the best known bound is v(k) <C fclog(fc) 
[H Chapter 12], which is derived from the usual Waring's problem. For small values 
of k, the best bounds for v(k) are derived from ideal solutions of the PTE problem. 
In fact, these are much better than those which derive from the usual Waring 
problem. See again [3, Chapter 12] for a full explanation of the connection between 
the two problems. 

In 1935, Wright [55] conjectured that ideal solutions should exist for all n. 
However, it does not appear that this conjecture is close to being resolved. For 
n = 2,3,4,5, complete parametric ideal solutions are known. For n = 6,7,8, only 
incomplete parametric solutions are known. See [3] Chapter 11] and [5] [U [9] for 
further details of these cases. For n = 10, infinite inequivalent families of solutions 
are known (albeit incomplete) [23) . 

For size 9, only two inequivalent solutions are known. These were found com- 
putationally by P. Borwein, Lisonek and Percival [4]. Until 2008, there were also 
only two inequivalent solutions known for size n = 12. They were both found com- 
putationally, by Kuosa, Myrignac and Shuwen [22] and Broadhurst [6]. However, 
in 2008, Choudhry and Wroblewski [11] found some infinite inequivalent families of 
solutions for n = 12 (again incomplete). Both infinite families of solutions for sizes 
10 and 12 arise from rational points on elliptic curves using a method of Letac's 
from 1934, which appears in [15] . 

For n = 11 and n > 13, no ideal solutions are known. 

Analytic methods are no closer to resolving Wright's conjecture. Along the same 
lines as the "Easier" Waring problem, define N(k) to be the least n such that the 
PTE problem of degree k has a solution of size n. Much work has been done on 
obtaining upper bounds for N(k), for example, see [13 [20] [26] and [3] Chapter 12]. 
The best upper bound is due to Mclzak, which is N(k) < \{k 2 — 3) when k is odd, 
and N(k) < \{k 2 —4) when k is even. Meanwhile, there is no lower bound on N(k) 
that would rule out ideal solutions. 

Although the PTE problem is traditionally looked at over Z, it may be viewed 
over any ring. Alpers and Tijdeman p] were the first to consider the PTE problem 
over a ring other than the integers. Their article discusses the PTE problem over 
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the ring Z x Z and shows that ideal solutions of size n in this case come from a 
particular kind of convex 2n-gons. Their article also gives an example of a solution 
to the PTE problem over the Gaussian integers, Z[i]. It further notes that there 
does not appear to be any other mention in the literature of the PTE problem in 
this setting. In [TO], Choudhry examines the PTE problem over the ring of 2 x 2 
integer matrices, A/2(Z). 

Based upon the work of Alpers and Tijdeman, because Z[i] contains Z, we might 
expect smaller ideal solutions to the PTE problem in this setting. Therefore, this 
article examines ideal solutions to the PTE problem over In particular, we 

view ideal solutions over Z as special cases of solutions over Z[i]. We also describe 
a computational search for ideal solutions of size n for Z[i] for n > 8. This search 
generalizes the methods of Borwein, Lisonek and Percival in [4]. They performed 
a computer search for ideal solutions of size n = 10 and n — 12, which took 
advantage of an alternative formulation of the problem to reduce the number of 
variables. Their search was further optimized by using the arithmetic properties of 
ideal solutions. 

All the results that are required for the method of Borwein et al. generalize 
sufficiently to Z[i]. We proceed by discussing some further background from the 
existing literature in Section [21 and then explaining the computational method to 
be used Section [3] In Section [H we will provide analogues of existing theorems 
in the literature for the PTE problem over the Gaussian integers, which allow the 
computational search to be optimized. Finally in Section [6j we describe the results 
of a computational search for ideal solutions for n — 10 and n = 12. 

For convenience, we state some results in greater generality than Z[i]. As a 
general notation, we refer to the PTE problem over the ring R as the i?-PTE problem. 
Throughout this article, let £ € C be an algebraic integer, and let O denote the ring 
of integers of the number field Q(£)- Note that it is easy to find O-PTE solutions, 
such as the example found by Goldbach given above. Hence, we proceed to discuss 
the PTE problem in this general setting. 



Solutions to the Z-PTE problem satisfy many relations. Most of them generalize 
to O in a completely trivial way, and can easily be proved using Newton's identities. 
We list a few of them. Suppose X = {x\, . . . , x n } and Y = {y\, . . . , y n } are subsets 
of O, and k £ N with k < n — 1. Then the following relations are equivalent: 



2. Background 



n 



n 



(2.1) 





i=l 



( = 1 



(2.2) 




(2.3) 




Note that for any c G C, we have z c 



e cln(z) 



Since 
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and wc merely need differentiation for the proof of (|2.2p <^=> ()2.3|) . we can use this 
fact formally. Similarly, since the terms in the sum X)"=i zXi ~ S™=i zVi are n °t> m 
general, polynomials, we consider the division in (12.31) to refer to the order of the 
zero at 1. These relations provide an alternative formulation for the PTE problem. 
Note that in particular, the relation (|2.1| <^=> (|2.2[) implies that when X = n -l Y, 



Y[{z-x l )-Y[(z~y i ) 

i=l i=l 

is a constant in O. This constant plays a significant role in the study of the PTE 
problem, which we discuss later in Sections |3] and SJ 

Given a solution to the O-PTE problem, we can generate an infinite family of solu- 
tions. That is, if {xi, . . . , x n } and {y%, . . . , y n } are subsets of O with {xi, . . . , x n } =k 
{yi,...,y n }, then 

(2.4) {Mx x +K,...,Mx n + K} = k {My, + K,..,,My n + K}, 

for any M, K € O. This fact leads us to give the following definition: 

Definition 2.1. Let Q(C) be a number field and O be its ring of intgers. Suppose 
Xi =fe Yi and X2 =k Y 2 . If there exists an affine transformation f(x) — Mx + K 
with M, K in Q(C) such that /(Xi) = X 2 and f{Yi) = Y 2 , then we say that 
Xi =fc Y\ and X 2 =kY 2 are equivalent. 

The following fact can be used as a criterion for PTE solutions to be equivalent. 
It will be useful later. 

Proposition 2.2. Suppose {x%, . . . ,x n } = n -i {j/i, • ■ • , Vn} and {x[, . . . , x' n } = n -i 
{y'l, . . . , y' n } are equivalent ideal PTE solutions via the transformation f(x) = Mx + 
K where M,K e Q(C). UUti( x ~ x i) ~ Uti( x ~ Vi) = C and Uti( x ~ O ~ 

nr=i(^ - y'd = c '> thm c> = cMn - 

Proof. If these solutions are equivalent, without loss generality, we may assume 
Mxi + K — x' { and My± + K = y^ for i = 1, . . . , n. Thus, we have 

n n 

l[{x - {MXi + K)) - l[(x - (My, + K)) = C , 
i=i i=i 

and since this holds for all values of a;, we may replace x by x + K to obtain 

n n 

Y[((x + K)- [Mxi + K)) - Y[((x + K) — [My, + K)) = C . 
i=i i=i 

Simplifying, dividing through by M n and then replacing x/M by x, we obtain 

J[(x - Xi ) - J[(x - Vi ) = — , 

i=l i=l 

proving the result. □ 

Let z denote the complex conjugate of z. It is clear that if {x\, . . . , x n } =k 
{Ui,---,Vn}, then {x x ,...,x n } = k {Fi, •••,¥„} also. 
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3. Searching for Ideal Solutions Computationally 

We might naively search for ideal solutions to the PTE problem over Z in the 
following way. Suppose our search space is Xi,yi <E [0,5] PI Z. We may assume 
xi = 0. Then select the remaining integers so that < x 2 < Xa < • • • < x n and 
1 < Vi < ■ ■ ■ < y n -i, and take y n = x\ + . . . + x n - (yi + ■ ■ ■ + y n -x)- Now check 
whether or not 

x\ + . . . + x k n = y\ + . . . + y\ 

for each k — 2, . . . ,n — 1. This method requires searching in 2n — 1 variables. 

However, Borwein et al. [4] improve on this significantly. Recall from (|2.2|) that 
if {xi, . . . , x n } = n -i {yi, . . . , y n } is an ideal PTE solution, then 

(z - xi)(z -x 2 )---(z- x n ) - (z - yi)(z -y2j"-{z- y n ) = C, 

for some constant C € Z. Rearranging this equation and then substituting z — yj 
for j = 1, . . . , n we obtain 

(3.1) (yj - xi) ■■■{yj- x n ) = C. 

For any k G {1, . . . , n}, equation (|3.ip can be rearranged to 

(3 - 2) ~ ^- fe+2) ' ' ' ~ ^ - ^-x^-l.-x^y 

Now define 

f(z) := —{z - x n - k+2 ) ■ ■ ■ (z - x n ). 

From (|32J|, we have f(yj) = (aj _ 2 . l) ...(^_- Cti _ fc+l) for j = 1, . . . , k. So if the vari- 
ables xi, . . . , Xn-k+i and yi, . . . ,yk are known, then we also have the ordered pairs 
(yj,f(yj)) for j = 1, . . . , k. We may determine f(z) uniquely by using Lagrange 
polynomials and the ordered pairs (yj, f(yj)) for j = 1, . . . , k (see, for example [UJ 
Chapter 5]). Thus, f(x) is a polynomial of degree k— 1, and solving /(z) = yields 
its roots, which are x n -k+2, ■ ■ ■ , x n . Repeating this process gives the remaining 

Vk+i, ■■■,y n - 

The method of Borwein et al. requires searching through only n + 1 variables, 
instead of 2n — 1. This method can clearly be generalized to any ring of integers 
O, and this is what we have implemented for O = Z[i]. 

3.1. Optimizing the Search. In order to explain how Borwein et al. further 
optimize the search over Z, we need a definition, and in order to state it, we now 
restrict ourselves to any O that is also a unique factorization domain (ufd). We 
maintain this restriction for the remainder of this article. 

Definition 3.1. Suppose O is a ufd. Suppose X = n -i Y is a O-pte solution with 

nr=i( z _ x i) _ n"=i(^ - y») = c n ,x,y Tnen iet 

C n := gcd{C n ,x,Y|X =n-i Y}. 
We say that C n is the constant associated with the O-pte problem of size n. 

Thus, C n keeps track of all the common factors that appear among the con- 
stants that come from the second formulation of the PTE problem in (|2.2|) . The 
requirement that O is a ufd is necessary for C n to be well-defined. 

We have the following Theorem, generalized from Proposition 3 in [4], 
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Theorem 3.2. Suppose O is a ufd. Let {xi, . . . , x n } = n -i {yi, • ■ • , Un} be subsets 
of O that are an ideal O-pte solution. Suppose that q 6 O is a prime such that 
q | C n . Then we can reorder the yi such that 

Xi = yi (mod q) for i = 1, . . . , n. 

The proof of Theorem 13.21 follows that of Proposition 3 from [3] , but we repeat 
it for completeness. 

Proof. Assume q G O is a prime dividing C n . Since O is an integral domain and 
q is prime, (q) is a prime ideal. Since prime ideals of rings of integers of number 
fields are also maximal (see, for example [23]), the quotient Of(q) is a field. Let 
¥ q denote this field. It follows that IlILi^ — x *) ~ TYi=i( z ~ Vi) ecma l s a constant 
times q, and so is zero in F g [z]. Hence, Iir=i( z — x i) ~ n™=i( z — Vi) m Since 
¥ q is a field, the polynomial ring ¥ q [z] is a unique factorization domain. Since each 
of the factors z — Xi and z — yi are irreducible, it follows that the sets {x%, . . . , x n } 
and {yi, . . . ,y n } are equal as subsets of ¥ q . That is, they are equal modulo q, as 
desired. □ 

Borwein et al. [4] use Theorem l3.2l to optimize the search for ideal solutions over 
Z. This can also be applied over O. Suppose q\,q2 are the two largest primes (in 
O) dividing C n . Assume x\ — 0, and pick the rest of the variables so that for 
i = 1, . . . , n 

Xi ee yi (mod qi) 

i 

- yi) ■ ^(xj - y 3 ) ee (mod q 2 ). 
j=i 

We assume the solutions Xi and yt pair modulo qi , and that each Xi pairs to the 
previous yi modulo q 2 , unless all the Xj and yj are already paired off modulo q 2 . 
Hence, every prime q that divides C n reduces the search space in each variable by a 
factor of N(q), where N(q) denotes the algebraic norm of q. Therefore, divisibility 
results, particularly large prime factors, for C n are very important for optimizing 
the search. 

4. Divisibility Results for C n 

There are a number of results in the literature concerning divisibility of C n for 
the Z-PTE problem. For example, about half of the article by Rees and Smyth [3T] 
is spent proving such results. Many of these results generalize immediately to O, 
which we state below without proof. In the case where the result is more of an 
analogy than a generalization, we provide a proof. 

The usual method of generalization is to view arithmetic modulo a prime power 
in Z as analogous to arithmetic in the appropriate finite field, which is then viewed 
as analogous to arithmetic modulo the algebraic norm of a prime in O. Fermat's 
Little Theorem corresponds with Lagrange's Theorem and so on. This method was 
used in the proof of Theorem l3.2l in the previous section. 

The next two results are generalizations of Proposition 2.3 and Proposition 3.1 
in [21], respectively. 

Theorem 4.1. Suppose O is a UFD. Let q E O be a prime with N{q) > 3. Then 
N(q)\C N{q) . 



THE PROUHET-TARRY-ESCOTT PROBLEM FOR GAUSSIAN INTEGERS 



7 



Theorem 4.2. Suppose O is a UFD. Let q £ O be a prime such that 

n — 2 

n + 3 < JV(g) < n + 3 + — — . 

Then q \ C n +i ■ 

Note that Rees and Smyth use a "Multiplicity Lemma" to prove this result in 
[2Tj . The proof of this lemma also generalizes appropriately to O, and so Theorem 
14.21 is remains valid. 

We now prove a general divisibility result of C n for powers of primes q. This 
result is based on the same techniques used in Proposition 2.4 of |21) . 

Proposition 4.3. Suppose O is a UFD, and q £ O is a prime. If q | C n , then 



<1 



\ «"<,) 1 



where \x\ denotes the smallest integer greater than x. 

Proof. Suppose X = {<zi,...,a„} and Y = {bi,...,b n } with X = n -i Y, and 
9 I Cn.A.Y- From Theorem 13.21 we can relabel the Oj and bj such that ai = bi 
(mod q) for i = 1, . . . ,n. Note that O has N(q) congruence classes modulo q, and 
so there is at least one congruence class with at least \n/N(q)~\ elements from the 
set {b±, ...,&„}. Relabel this set so that b\, . . . , b |- n i are in the same congruence 
class modulo q. From equation (|2.4p . we can shift the at and bi by — b\, giving 
C n ,x,Y = ai«2 ■ • • a n . Then 

ai = bi = (mod q) 
a-i = bi = (mod q) 



fir n i = fer a i = (mod o). 

JV(g) I I N(q) I V " 

Thus, g^iv^yl | C n ,x,Y, and since X and T were arbitrary, we have proved the 
result. □ 

Note that we can only apply Proposition ^. 31 when we already have from another 
source that p \ C n . 

We now prove a specific result for the divisibility of C5 for powers of primes q G 0, 
with N(q) — 2. This result is based on the same techniques used in Proposition 2.5 
of HQ. 

Proposition 4.4. Suppose q £ O is prime with N(q) = 2. Then q A \ C5. 

Proof. Suppose X — {a\, . . . , a§} and Y — {61, . . . , 65} with X =4 Y . As in the 
proof of Proposition 14.31 we can relabel the ai and bj such that ai = bi (mod q) 
for i = 1, . . . , 5, and so that b%, . . . , 63 are in the same congruence class modulo q. 
Again as above, we can shift the 04 and &i by —61, giving C^,x,y = 0102030405. 
Assume that g 4 j C$ y x.Y- Since we know that g 3 | C^^xy however, we can assume 
that ai = 02 = 03 = g (mod g 2 ) and 04 = 05 = 1 (mod g). As usual, we have 

(4.1) (z-a 1 )(z-a 2 )(z~a 3 )(z~a 4: )(z-a 5 )-z(z-b2)(z-b 3 )(z-b4 : )(z-b 5 ) = C 5iX ,y- 

Substituting z = a% into (|4.ip gives 

-ai(ai - o 2 )(ai - b 3 ) (ax - &4)(ax - 65) = C 5< x,y- 
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Since a\,a\ — 62, a>i — 63 are all equivalent to modulo g, while a\ — 64, a\ — 65 are 
both equivalent to 1 modulo g and their product is not divisible by g 4 , we must 
have ai = a\ — 62 = a>i — 63 = g (mod q 2 ). Since a\ = q (mod q 2 ) already, this 
means that 62 = 63 = (mod q 2 ). 

We now substitute z = 04 into (14. ip giving 

-a 4 (a 4 - 6 2 )(a 4 - 6 3 )(a 4 - 6 4 )(a 4 - 65) = C 5 ,x,Y- 

Since 04,04 — 62,04 — 63 are all equivalent to 1 modulo q, while 04 — 64,04 — 65 
are both equivalent to modulo q, we can assume, without loss of generality, that 
04 — 65 = q (mod q 2 ) and 04 — 64 = g 2 (mod g 3 ), i.e., 04 — 64 = (mod q 2 ). 
Finally, substituting x — 65 into (|4.1I) gives 

(65 - ai)(6 5 - a 2 )(6 5 - 03X65 - 04) (65 - a 5 ) = C 5j x,y- 

Only &5 — 04 and 65 — 05 are equivalent to modulo q. However, we already have 
that CI4 — 65 = q (mod q 2 ), and so we must have 05 — 65 = (mod g 2 ). However, 



Not all results from the literature concerning C n generalize to Z[i] or O, and 
some must be addressed specifically depending on the ring of integers involved. 
Those relevant to our computer search for ideal solutions over Z[i] are discussed 
next. 



An important divisibility result for C n over Z is that n! | C n +i (see Proposition 
2.1 in [5T|, originally due to H. Kleiman in [IB])- This fact demonstrates that C n 
is highly composite and will contain some large prime factors. The proof that Rees 
and Smyth provide of Proposition 2.1 in [3T] uses the obvious fact that if t € Z then 
t(t + l)(t + 2) . . . (t + n) = (mod (n + 1)!). However, this depends on t being an 
integer. Unfortunately, this fact does not fully generalize to Nevertheless, we 
are able to state an analogous lemma below. For completeness, we prove this result 
in greater generality than necessary, that is, for the ring of integers of an arbitrary 
quadratic number field. 

We first recall some facts concerning quadratic number fields from Chapter 5 of 
[7]. Let K = Q(vd) be a quadratic field with d 7^ 1 squarefree and let D = d(K) 
denote the discriminant of K, and let O be its ring of integers. We also assume 
that O is a UFD, but note that this hypothesis is not required for Propositions 15.11 
and 15.21 or Lemma T5. 31 We have the following results: 

Proposition 5.1. (i) If d = 1 (mod 4), then {1, 1+ 2 } * s an integral basis 
for O and D = d. 

(ii) If d = 2 or 3 (mod 4), then {1, V^} is an integral basis for O and D = 4d. 
Thus, we may write O = Z[w], where uj = ■ 



Proposition 5.2. Let p be a prime and ( |J be the Legendre symbol. Then the 
decomposition of prime ideals of 1 in O is as follows: 



we have 



= ai + 02 + 03 + a 4 + 05 — (6 2 + 03 + 64 + 05) 
= 9 + 5 + + 64 + 65 — — — 64-65 = (mod g 2 ), 
which is a contradiction, proving the proposition. 



□ 



5. Divisibility Results for C n over Z[i] 
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(i) If p | D, i.e., if \ = 0, then p is ramified, and we have pO — p 2 , where 
p = pO + ujO, except when p = 2 and D = 12 (mod 16). In this case 
p = pO + (1 + u)0. 

(ii) If (j^J = — 1) then p is inert, and hence p = pO is a prime ideal. 

(iii) If (|r) = ^' fien P is sp/ii, and we have pO = pip2,. where pi = + 

(w — -^r^) O and p2 = pO + (ui — and c is any solution to the con- 

gruence c 2 = D (mod 4p). 

Lemma 5.3. Let p G Z &e a prime that is either ramified or split in O, i.e., is of 
type (i) or (iii) from Proposition \5.2\ Let s G N. Then 

where p is type (i) and pO = p 2 , 
t(t + l)(t + 2)(t + 3) . . . (t + sp- 1) G { pf where p is type (iii) and pO — pip 2 , 

where p is type (iii) and pO = pip2- 

Proof. Define a map <fr : O —> M 2 by <j>(a + bui) = (a, b), where a, b G Z. Because 
{1, ui} is an integral basis for O, it is clear that <j> is well defined. We now examine 
the image of ramified and split ideals pO under <f>. 

First note that uj satisfies the equation ui 2 = D ~P — h Duj. 

Suppose p is ramified. Then from Proposition 15.21 we have pO = p 2 , where 
p = pO + uiO, excluding the case that p = 2 and D = 12 (mod 16). Thus, an 
arbitrary element q G p looks like 

q = p{a + bui) + u)(e + fuj) 
— ap + (bp + e)uj + fui 2 

(D-D 2 

= ap+ (bp + e)u + f I h Duj 

fD-D 2 \ , 
= ap + f( j +(bp + e + Df)uj, 

where a, b,e, f G Z. Thus, we have 

4>(q) = ( ffl P + /( £> ^ £>2 ) >bp + e + Df 

In the case that p = 2 and D = 12 (mod 16), we have p = pO + (l+uj)0. Thus, 
an arbitrary element q G p looks like 

g = p(a + bui) + (1 + w)(e + /w) 
= + e + (6p + e + f)ui + fui 2 

(D-D 2 

= ap + e + (bp + e + f)u + f I h Duj 

(D-D 2 \ 
= ap + e + f ( J + (bp + e + (D + 

where a, 6, e, / G Z. Thus, we have 

0(«)= (ap + e + f( — A ° 2 ) ,bp + e + (D + l)f 
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Alternatively, suppose p is split. Then from Proposition l5.21 we have pO — pip 2 , 

O and p 2 = pO + (u - ^ 



where pi = pO + (w — ^t 2 ) and p 2 = pO + (u) — -^r 2 ) and c is any solution to 



the congruence c 2 = D (mod 4p). 

Thus, an arbitrary element of q € pi (resp. p 2 ) looks like 

, , / Z3±c\ , 
9 = P(a + buj) + I u — I (e + /w) 

= ap + 6pw + ea; + /cj 2 + / I + Doj J - I — ^— J e + I — ^— ) fu 

_, + ,(^)_(^). + (* +e+/D+ (^),)„/ 

where a, b,e, f £ Z. Thus, we have 

«„ = („ + f (^) - (£±£) e , („ + . + , D + (£±£ ) , 

Now let i 6 and suppose i = u + wcj so that <p(t) = (u, v). Note that D ± c 

is always even, and if we pick / so that / ( °~® j is an integer, then 0(g) e Z 2 

in all three of the above cases. Further, in each case, we may solve the equations 
bp + e + Df = v, bp + e + (D + l)f = v and bp + e + fD + (£±£) / = v for 6, e, /. 
Thus, in each case, it follows that the set 

{t + jp,t + jp + l,t + jp + 2,t + jp + 3,...,t + jp + (p-l)} 

contains an element that belongs to p or p L or p 2 respectively, for j = 0, . . . , s — 1. 
Thus, the set {t, t + 1, t + 2, . . . , t + sp — 1} contains s elements that belong to p or 
pi or p 2 respectively, proving the lemma. □ 

Remark 5.4. Note that if p is inert, i.e. of type (ii), the expression t(t+ l)(i + 2)(i + 
3) ...(£+ n) need not belong to pO. For example, when K = Q(i) and O = Z[i], 
p = 3 and t = i, note that none of i, 1 + i, 2 + i, . . . , n + i contain a factor of 3. 

Using the above characterization of primes in a quadratic number field, we have 
the following result for the O-PTE problem analogous to nl | C n+ i. 

Theorem 5.5. Let C n +i be the constant associated with the O-PTE problem. Sup- 
pose p is either (i) ramified or (Hi) split and let 

where p is type (i) and pO = p 2 , 
p = ^ pi where p is type (Hi) and pO = pip 2 , 
where p is type (Hi) and pO = pip 2 

Let s = \_(n + l)/pj and let £ be the highest power such that n + 1 6 p £ . Then 

Cr- t-max(s— £,0) 
n+1 £ P ^ • 

We digress before proving Theorem 15.51 As stated earlier, many results on the 
O-PTE problem involve Newton's identities and symmetric polynomials, including 
the proof of (|2.ip <^=> (|2.2[) . We need them for the proof of some results below, so 
although they are well known and easily found in the literature (for example see 
[24]), we repeat them here. 

Let ntN. Let si, . . . ,s n be variables. Then for all integers k > 1, we define 

p k (s!, s n ) := s* + s 2 + . . . + s*, 
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the fcth power sum in n variables. Similarly, for k > 0, we define 



e (si, . . 


■i s n) 


= 1 


ei(si, . . 


■i s n) 


= Si + S 2 + ■ 


e 2 (si, . . 


■i s n) 


= SiSj 






i<j 


e„(si, . . 


■i s n) 


= sis 2 ■■ ■ s n 


e fe (si, . . 


■ ) $n) 


= 0,Vfc > n, 



to be the elementary symmetric polynomials in n variables. Then we have the 
result known as Newton's identities: 

k 

(5.1) ke k (si,...,s n ) = y^-iy^efc^si, . . . , s n )Pi(si, ■ • - , s n ), 

i=l 

for all k > 1. Note that this can be rearranged to 
(5.2) 

fc-i 

Pfc(si, . . . , s„) = (-l) fe_1 fce fc (si, . . . , s„) + ^(-l)' £_l e fe _ i (si, . . . , s„)p i (si, . . . , a„), 

i=l 

for k > 2. Another fact is the identity 

n n 

(5.3) J]^ - s ^ = I^- 1 )*^*!, ■ ■ ■ , s n )t n - k . 

i=l fc=0 

Thus, the coefficients of a polynomial are elementary symmetric polynomials of its 
roots, and because of (|5.ip . they depend on the power sums Pi(si, . . . , s n ). 
We are now able to proceed with the proof of the Theorem. 

Proof of Theorem ] 5. 51 We closely emulate the proof of Proposition 2.1 in [2T|. Sup- 
pose X = {xi, . . . ,x n +i} and Y = {yi, . . . ,y n +i} ar e subsets of O, and X = n Y. 
Then we have 

(z - xi) ■ ■ ■ (z - x n+ i) - (z - yi) ■ ■ ■ (z - y n +i) = C n+ i >x ,Y 

where C n +i,x,Y = (-l) n+1 (a;ia;2 • • • x n +i - y\y% • ■ ■ Vn+i)- Then from the identity 
(|5.3|) . we have 

n+l 

{z- Xl )---(z- x n+ i) = ^2(-l) k e k (x u . . . ,x n+1 )z n+1 ~ k 

k=0 
n+l 

{z- yi )---(z- y n+1 ) = ^(-l) fe e fc ( yi , . . . , y n+ i)z n+1 - k . 

k=0 

From the identity (|5.2| . it follows that 

(5.4) pk+i(xi, ■ ■ ■ ,x n +i) = (-1) (k + l)e fc+ i(xi, . . .,x n +i) 

k 



y^(-l) fe+1 l e k +i-i(x!, . . .,x n+ t)pi(xi, x n+ i), 



i=i 
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and 

(5.5) Pk+i(yi, ■ ■ .,y n +i) = (-l) fe (fc + l)efc+i(yi, . . • ,jm+i) 

k 

+ ^(-l) fe+1_, e i;+ i_ i (y 1 , . . . ,y n +i)Pi(yi, ■ ■ • ,2M+i)- 

i=l 

By hypothesis, we ha.vep k (x 1 , . . . , x n+1 ) = p k (y 1 ,..., y n+1 ) and e k (xi,..., x n+1 ) ■ 
efc(yi, . . . , Vn+i) f° r 1 < k < n, and so subtracting (|5.5|) from (|5.4p it follows that 

p n+ i(xi,...,a; n+ i) -p„+i(yi, . . . ,2/„ + i) = (-l) n (n + l)e n+ i(xi, . . . , x n+ i) 

(5.6) - (-l) n (n + l)e„+i(yi, . . . ,x y+l ). 

Now noting that C„+i,x,y = e n+ \(x\,. . .,x n +\ - e n+ i)(yi, . . . ,y„+i), rearranging 
(|5.6|) we get 

(5.7) p n+ i(ari, . . . ,x„+i) + (-l) n (n + l)Cn+i,x,Y = p„+i(j/i, • . .,y n+ i). 

Since s = [(n + l)/p\ , it follows that sp < n + 2, and so from the above lemma 
we have 

n+l 

(5.8) ^(-l) fc e fe (0, -1, . . . , -n)t n+1 - k = t(t + l)(t + 2)(t + 3) . . . (t + n) € p s . 

fc=0 

Substituting t = x\,X2, ■ ■ ■ ,x n +i into (|5 .8[) and summing, and doing the same 
for t = j/i, t/2, • • • , J/n+i, we get 



n+l n+l 
= 1 fe=0 



n+l n+l 



( 5 - 9 ) E E^ 1 )^^ - 1 - • • • > e p s 

and 

( 5 - 10 ) E Et- 1 )^ * - 1 ' • • - -™)^ +1 ~ fe e p s . 

Subtracting (|5.10l) from (|5.9[) and applying (|5.7p . we get 

(n + l)C n+ltX ,Y G p s . 

Since n+l ep f and n + 1 ^ p £+1 , we have C n+ i, x ,y G p max ( s -^ ) ; an d because X 
and Y were arbitrary solutions to the O-pte problem, we have C n +i € pmax(s-«,o)^ 
proving the theorem. □ 

The above divisibility results give lower bounds for C n for the PTE problem over 
Z[z]; these are stated in Table [TJ 

When (*) appears in Table [TJ this means the upper bounds for C n come from 
the upper bounds for C n for the PTE problem over Z (which we have included for 
comparison in Table [5]) . In the next section we explain the upper bounds new 
to the Z[i]-PTE problem. These have been determined by searching for solutions 
computationally. 

Table [2j which largely comes from [4 , 5 but is updated with the new solutions 
from [6] and [UJ, shows that the constant for the Z-PTE problem has many more 
factors than that for the Z[i]- PTE problem. This demonstrates that the Gaussian 
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Table I. Divisibility Results for the Z[z]-pte Problem 



n 


lower bound for C„ 


upper bound for C n 


2 


1 


1 


3 


(i + a 


(i + a 


4 


1 


1 


5 


(l + i) 4 (2 + i)(2-i) 


(l + z)°(2 + z)(2-z) 


6 


(1 + i) 3 (2 + i)(2 - i) 




7 


(l + i) 4 (2 + i)(2-i) -3 


(l + i) (2-i)^(2 + i)^3 


8 




Cl + i) s C2 - iY(2 + i) 2 C3 + 2i)(3 - 2i) 


9 


(l + i)°(2 + i)(2-z) 
•3 2 • (3 + 2i)(3 - 2i) 


(1 + i) LS {2 - if{2 + if ■ 3 4 • I' ■ 11 • (3 + 2i) 
(-3 + 2z)(4 + z')(4 - i) ■ 23 • (5 + 2i)(5 - 2i) (*) 


10 


(l + i) 5 (2 + i)(2-i) 
(3 + 2i)(3-2i) 


(1 + z) i3 (2 - i) a (2 + i) 2 • 3 2 • (3 + 2i)(-3 + 2i) 
(4 + z)(4-z) 


11 


(l + i)"(2 + i)^(2-i)^ 


none known 


12 


(l + i)°(2 + i)^(2~i)^ 


(1 + i)" 4 (2 - i) J (2 + if ■ 3 a • 7* ■ 11* • (3 + 2if 
(-3 + 2i) 2 (4 + z)(4 - z) • 19 • 23 • (5 + 2i) 
(5 - 2z) • 31 (*) 


13 


(l + t)''(2 + t) : '(2-i) il 

(3 + 2»)(3-2t)(4 + t)(4-t) 


none known 


14 


(l + i)'(2 + i) ii (2-i)^ 

(3 + 2i)(3-2i)(4 + z)(4-i) 


none known 


15 


(l + i) 8 (2 + z)(2-i) 
(3 + 2i)(3-2i) 


none known 



Table 2. Divisibility Results for the Z-PTE Problem 



n 


Lower bound for C n / (n — 1)! 


Upper bound for C n /(n — 1)1 


2 


1 


1 


3 


2 


2 


4 


2 • 3 


2 • 3 


5 


2-3-5 


2-3-5 


6 


2 2 • 3 • 5 


2 3 • 3 • 5 


7 


3-5- 7- 11 


2 2 • 3 -5 -7- 11 -19 


8 


3-5- 7- 11 


2 4 - 3 -5 -7- 11 -13 


9 


3-5- 7- 11 


2 2 -3 2 -5-7- 11 -13- 17-23-29 


10 


5-7-13 


2 4 -3 2 -5-7-ll-13-17-23-37-53-61 
•79 -83 -103 -107 -109 -113 -191 


11 


5 -7- 11 -13 -17 


none known 


12 


5-7-11 


2 4 • 3 5 • 5 • 7 • 11 • 13 2 • 17 • 19 • 23 • 29 • 31 



integers are a much less restrictive setting for the PTE-problem than the ordinary 
integers. 

6. Computer Search for Solutions 

We may restrict the O-pte problem to a symmetric version. This is helpful 
because there are fewer variables, but at the same time, some ideal solutions 
may be missed. For odd n, this means finding solutions xi , . . . , Xn e O with 
{xi, . . . , x n } = n -i {—xi, . . . , — x n }- Since xf k = (—Xi) 2k for all k e N, this means 
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we only need to consider solutions to E<=i x t ~ f° r e = 1, 3, . . . , n — 2. For exam- 
ple, {3+3i,3+4i,3+5i,-2-8i,-7-4i} = 4 {— 3 — 3i, -3-4i, -3-5i, 2 + 8i, 7 + 4i} 
is an ideal symmetric solution of size 5. 

Similarly, for even n, this means finding solutions x%, . . . , a: n /2, yi, . . . , y ra /2 G O 
such that {xi,...,a; n/ 2,-a;i,...,-x„ /2 } = n _i {y x , . . . ,y n/2 , -y 1 , . . . , -y n/2 } . As 
above, since (— Xi) 2e+1 = —x 2e+1 for all k € N, we only need to consider solutions to 

E?ii = E?ii 2/f for e = 2,4, . . . , n- 2. For example, {±1, ±(4 + i), ±(3 + *)} =5 
{±(7i), ±(7 + 4i), ±(7 — 3«)} is an ideal symmetric solution of size 6. 

Thus, the symmetric case of PTE problem involves half as many variables as the 
usual case. Some results concerning this case are discussed in 0] [5J [§] ■ 

In [1] , Borwein et al. describe an algorithm for finding odd and even symmetric 
solutions to the Z-PTE problem. We have adapted this algorithm for finding ordi- 
nary solutions as well as odd and even symmetric solutions to the Z[i]-pte problem. 
As the ideas behind the algorithms are not any different from the original, one may 
see [4] for an explanation. This was implemented first in Maple and then in C++, 
using the Class Library for Numbers. 

The computer search was implemented to try to find solutions with real and 
imaginary parts between and 30 for sizes 10 and 12. The above method is trivially 
parallelizable, so each search range was divided up into intervals, which were then 
submitted to a cluster of machines. 

The following symmetric solutions of size 10 were found: 

{±(9 + i), ± (4 + 8i), ±(8 + U), ±(3 - 3i), ±(1 - 9i)} = g 
{ ± (5 + 7i), ±8, ±(9 + 3i), ±8i, ±(1 - 7i)} 

which has constant 

-(1 + i) 22 (2 + if {2 - i) 2 3 2 (3 + 2i) 2 (3 - 2i)(4 + i){A - i)(5 + 2i), 

and also 

{±(8 + 3i), ± (9 + 4i), ±(11 + 2t), ±(1 - 7<), ±(5 + 7*)} = 9 

{ ± (7 + 7t),±(ll + H),±(ll + 4i),±(l + 6i),±5i} 

which has constant 

i(l + i) 13 (2 + i) 2 {2 - i) 2 3 2 (3 + 2i) 2 (3 - 2i)(4 + i)(4 - i)(5 + 2i)(5 - 2i)(5 + 4«). 

Additionally, by the remark at the end of Section 2, the complex conjugates of 
these solutions are also ideal PTE solutions of size 10. First note that none of these 
solutions lies on a line in the complex plane. Thus they cannot be equivalent to a 
Z-PTE solution. By examining their constants and applying Proposition 12. 2[ they 
cannot be equivalent to each other either. 

Further note that all the Gaussian integers in the first solution have norm < 90, 
while in the second they all have norm < 147. This contrasts with the ordinary 
integer case where from [3] there is no size 10 solution with height less than 313, 
and in fact, there are only two inequivalent solutions with height less than 1500. 
This results corresponds to the intuition that Gaussian integer solutions should be 
"easier" to find. 

We now explain the second column of Table 1 above, which lists the upper bounds 
for the divisibility of C n . 

For n — 2 and n — 3, the upper bound comes from Table 2. 

For n = 4, the upper bound comes from the solution {0, 0, 0, 0} =3 {1, —1, i, —i}. 
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For n = 5, the solutions 

{0, -hi, -3 - 4i, 1 + 3i, 1 + 3i} = 4 {-5 - 5i, 5, -4 + 3i, 1 - 7i, 2 + 6i} 

and 

{0, 2 - 4i, 3 - i, -6 - 3i, -4 - 7i} = 4 {-5 - 5i, -4 - 2i, 4 - 3i, -2 - 6i, 2 + i} 

have constants — (1 + «) 5 (2 — i) 2 (2 + z) 7 and i(l + i) 6 (2-i)(2 + if respectively. The 
upper bound comes from taking the gcd of these constants, along with the constant 
associated to the complex conjugate of the second solution. 
For n = 6, the solution 

{0, -hi, 2-4i, -4-2i, -6+2i, -4+3», — 5— 5i} = 5 {-h+hi, l+3i, -8+i, l-7i, 4-3i} 

has constant — (l + i) 4 (2 — i) 2 (2 + i) 8 (— 3 + 2i). The upper bound comes from taking 
the gcd of this constant and the constant associated to the complex conjugate of 
this solution. 

For n — 7, the solution 

{3 + i, 2 + Ai, - 3 - 4i, 2 - 3i, -5 + 2i, -5 + 3i, 6 - 3i} = 6 

{ - 3 - i, -2 - 4i, 3 + 4«, -2 + 3i, 5 - 2i, 5 - 3i, -6 + 3i}, 

has constant (-i)(l + i) 6 (2 - z) 3 (2 + z) 2 3(3 + 2i)(4 + i)(5 - 2i). The upper bound 
comes from taking the gcd of this constant and the constant associated to the 
complex conjugate of this solution. 
For n = 8, the symmetric solution 

{±(2 + 2i),±3, ±3i, ±(2 - 2i)} = 7 {±2, ±2i, ±i, ±1} 

has constant (1 + if (2 - i) 2 (2 + i) 2 (3 + 2i)(— 3 + 2i). 

For n = 10, the upper bound is obtained taking the gcd of the constants from 
the solutions listed above, as well as their complex conjugates. 

For n = 9 and n = 12, the upper bound is obtained from factoring the bounds 
listed in Table 2. 

Note that for n in the range 4 < n < 8, many other Z[i]-PTE solutions are known, 
but they give no further information about the divisibility of C n . In these cases, 
we have not been able to prove if these upperbounds are true in general. 

Unfortunately, no symmetric solutions of size 12 have been found. Considering 
that there is a symmetric solution of size 12 of height 151 in the integer case, 
the usual intuition implies that a Gaussian integer solution would not be much 
larger than the search range. However, the search for size 12 ideal solutions took 
approximately 2 weeks on a cluster of 16 machines each with four lGhz. processors. 
Considering the magnitude of the solutions found in the integer case, this method 
does not seem likely to produce them in the Gaussian integer case. 

7. Further Work 

There are some natural directions for further work in this area. Clearly, the 
search range could be extended with the aim of finding more ideal solutions over 
the Gaussian integers. Additionally, since Theorem 15.51 (as well as the results in 
Section |4| holds for any ring of integers of a quadratic number field that is a unique 
factorization domain, a computer search as described in the previous section is 
certainly possible in any such ring. 
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However, our current implementation does not readily generalize to any ring of 
integers. Further, we believe that any implementation of a computer search in a 
different ring of integers would be significantly computationally slower than that for 
Z[t]. This is because for most mathematical software, determining whether or not 
a complex number is a Gaussian integer is naturally much easier than determining 
whether or not it is, say, an element of Z[e 27 ™]. 

We are currently working on another computational method which would gen- 
eralize the work of Smyth in [53] and Choudhry and Wroblewski in [TT]. Both of 
these articles relate the Z-PTE problem to elliptic curves. 
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